First-principles calculations of the magnetic properties of (Cd,Mn)Te nanocrystals 
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We investigate the electronic and magnetic properties of Mn-doped CdTe nanocrystals (NCs) 
with ^ 2 nm in diameter which can be experimentally synthesized with Mn atoms inside. Using the 
density functional theory, we consider two doping cases: NCs containing one or two Mn impurities. 
Although the Mn d peaks carry five up electrons in the dot, the local magnetic moment on the Mn 
site is 4.65 /is. It is smaller than 5 /is because of the sp-d hybridization between the localized 3d 
electrons of the Mn atoms and the s- and p-type valence states of the host compound. The sp-d 
hybridization induces small magnetic moments on the Mn nearest-neighbor Te sites, antiparallel to 
Mn moment affecting the p-type valence states of the undoped dot, as usual for a kinetic mediated 
exchange magnetic coupling. Furthermore, we calculate the parameters standing for the sp-d ex- 
change interactions: Conduction A^o^ and valence No (3 are close to the experimental bulk values 
when the Mn impurities occupy bulk-like NCs' central positions, and they tend to zero close to the 
surface. This behavior is further explained by an analysis of valence band-edge states showing that 
symmetry breaking splits the states and in consequence reduces the exchange. For two Mn atoms 
in several positions, the valence edge states show a further departure from an interpretation based 
in a perturbative treatment. We also calculate the d-d exchange interactions \J^^\ between Mn 
spins. The largest \J^^\ value is also for Mn atoms on bulk-like central sites; in comparison with 
the experimental d-d exchange constant in bulk Cdo.95Mno.05Te it is four times smaller. 



I. INTRODUCTION 

The recent progress in chemical synthesis, computa- 
tional capabilities and scanning-probe techniques has 
permitted a detailed understanding of semiconductor 
nanocrystals (NCs) - also known as nanoparticles, clus- 
ters, crystallites or quantum dots (QDs). Their prop- 
erties which depend on size^^ and shape^ have overim- 
posed effects due to quantum confinement. For instance, 
when the QD radius is smaller than the Bohr radius 
of exciton, quantum-confinement effects appear, such as 
the blue shift of gaps and the discretization of energy 
spectra^^. The gap properties can also be tuned inten- 
tionally with doping^—. Recently, much effort has fo- 
cused on II- VI semiconductor NCs doped with magnetic 
impurities such as Mn, which is the topic of this work. 
The Mn doping is motivated by diluted magnetic semi- 
conductors (DMSs). These compounds are bulk semicon- 
ductors doped with transition-metal impurities such as 
Cr, Mn, Fe and Co at low concentrations. DMSs are cur- 
rent research materials in spintronics, integrated in novel 
magnetoelectronic devices such as spin-LEDs^^. Their 
remarkable magnetic and magneto-optical properties re- 
sult from the strong sp-d exchange interactions between 
band carriers and Mn ions^^. These interactions yield gi- 
ant band-edge splittings at low temperature (about 100 
meV)i^. In particular, little is known theoretically about 
exchange in NCs made of II- VI DMSs doped with Mn. 

DMS NCs of type II- VI doped with Mn have been suc- 
cessfully synthesized and characterized during the last 
fifteen years. These works found at zero field, or low 
fields, contradictory results for the Zeeman splitting. In 



comparison with bulks, several authors have shown a re- 
duction of the excitonic Zeeman splitting in Mn-doped 
CdSl^ CdSeil and CdTe^^^^^^^^ NCs. Other works re- 
vealed no enhancement of the excitonic Zeeman splitting 
in Co-doped ZnO, ZnSe and CdS QDs^^ in comparison 
with their bulk values. Magnetic circular dichroism and 
optical experiments with Mn-doped ZnSe NCs revealed a 
value for the same splitting of 28 meV, much larger than 
in bulk ZnSe:Mni^ 

The Zeeman splittings are due on one hand to the 
confinement-enhanced sp-d hybridization between the 
occupied Mn 3d orbitals and the sp-type valence states of 
the host compound; and on the other hand to the crystal 
field experienced by Mn impurities, which is significantly 
different nearer the NC surface than in bulk^^. Experi- 
mentally, several authors found that Mn atoms are em- 
bedded in ZnSe NCs^^, and even a single Mn impurity 
is inside CdTe nanocrystals^^. Theoretically, enhanced 
splittings are obtained in ZnSe:Mn NCs within the effec- 
tive mass approximation— with a single Mn sitting at the 
center. However, there is little theoretical information 
about exchange interactions and splittings for CdTe:Mn 
NCs, and even less with Mn off-center. We also note that 
the DMS NCs can hold inside several Mn impurities. 

In the present work we study Mn-doped CdTe NCs of 
spherical shape with the density functional theory. In 
section [ni we give a brief account of the theorical frame- 
work and numerical details. The host compound is a 
well-known wide-gap semiconductor of type II- VI, and 
manganese is a widely-used dopant, known for activating 
photo- and electroluminescence, and also for contributing 
to efficient luminescence centers {3d electrons) Sev- 
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eral authors have previously calculated the properties of 
bulk II- VI CdTe doped with Mn^^^^. We consider two 
doping situations: NCs including one and two Mn impu- 
rities, which substitute for one or two Cd atoms in the 
zinc blende lattice, respectively. In section ITTTl the results 
concerning NCs with a single Mn impurity are given. 
The results for NC with two impurities are described 
in section IIVI We calculate the total ground-state en- 
ergy of the QD, for several positions of impurities within 
the crystal and the magnetic state. Moreover, we ob- 
tain the sp-d and d-d exchange constants^- as a function 
of the Mn locations: N^a and NqP stand for the ex- 
change interactions between the Mn local moments and 
the s- and p-type band-edge states of the host CdTe; J^^ 
parametrizes the exchange interaction between two Mn 
spins. Finally, we sum up the main findings and con- 
cluding remarks in section |Vl The system described here 
may provide further understanding of solid-state qubits, 
since it permits to detect and manipulate a single spin^^. 
In addition it could show magnetic and magneto-optical 
properties, such as fast recombination and high lumines- 
cence efficiency^. 



II. COMPUTATIONAL DETAILS 

The many-body problem for the electrons around the 
nuclei is solved based on the density functional theory 
(DFT), using the Kohn-Sham equations. The valence 
electrons move in the external potential created by the 
nuclei and the core electrons. The electronic states of the 
studied NCs are obtained from the projector augmented- 
wave method as implemented in the VASP code^^^^. 
To account for the sp-d hybridization^^ the ten 3c? spin- 
orbitals of Mn impurities are included in our calculations. 
Within DFT, our approach for the exchange-correlation 
is defined with the generalized-gradient approximation 
(GGA) of Perdew, Burke and Ernzerhof^^. However, 
the interactions among the 3 c? electrons of Mn impurities 
in the GGA approximation are only partially described, 
since they are strongly localized. Thus, we use the so- 
called GGA+ U scheme and introduce in the calculations 
two common correction parameters, U and J^^. For Mn 
atoms in CdTe bulk, = 6.2 eV and J = 0.86 eV. They 
are considered to correct the Coulomb ( U) and exchange 
(J) interactions among the 3(i electrons at Mn site^i^. 
We have checked that with the U values of Refs. |32| and 
[33I for Mn, the density of states are very similar, espe- 
cially the d peak corresponding to Mn is nearly at the 
same energy. This peak is below the p-type part of the 
CdTe valence band. We have chosen the recent U value 
in Ref. ^33 because it calculated MnTe bulk where Mn 
is surrounded by Te. Additionally, tests on bulk MnTe 
with the same lattice constant of Ref. 33| show nearly the 
same band structures when using LSDA or GGA and the 
same U value. 

The input parameters for the VASP calculations are 
determined in a preliminary work for two model bulk 



systems, CdTe and F-MnTe (zinc blende). For several 
lattice constants we fitted the ground-state energies to 
the Murnaghan's equation of state, which depends on 
the unit cell volume. The equilibrium lattice constants 
of bulks CdTe and F-MnTe are acdTe = 6.63 A and 
ctF-MnTe = 6.38 A. The cut-off energy in the plane 
wave basis set is 350 eV; this is the shared cut-off to 
converge the ground-state energies of both CdTe and F- 
MnTe bulks within meV. 

The studied NCs are spherical about 17 A in diameter 
centered on a cation site (Cd). The crystal structure is 
zinc blende, which has tetrahedral (Td) symmetry around 
the atoms. These QDs are passivated since otherwise the 
surface dangling bonds would introduce surface states in 
the near-gap spectrum. Organic ligands are commonly 
used to passivate NCs, but we want to center on the in- 
trinsic properties of NCs^^. For the sake of simplicity 
we resort to the simplest passivation agent: a pseudo- 
hydrogen atom (H*). The fictitious H* is characterized 
by fractional electronic charge and by a proton with the 
same charge but positive. Every Cd (5^^) dangling bond 
includes 2/4 unpaired electrons, so it is bound to a pseu- 
dohydrogen atom with a fictitious charge of 6/4 e. Sim- 
ilarly, every Te dangling bond (55^ 5p^) is bound to a 
pseudohydrogen atom with an electronic charge of 2/4 
e. When the pseudohydrogen atoms H* are included the 
NCs have 107 atoms in total: 19 Cd, 28 Te and 60 H*. 
We dope them with one or two Mn impurities depending 
on the studied case. 

We use the supercell approximation to deal with a sin- 
gle nanocrystal. We place the NC at the center of a 
large cubic unit cell and take the F point. Then we have 
converged the ground-state energies of two small model 
NCs ( CdTe4lI*i2 and MnTe4H*i2) versus the superceh 
size within me Vs. We consider that this separation be- 
tween walls is valid for larger NCs. The input distances 
between nearest-neighbor atoms are taken from those in 
Cd-Te bulk. We have dcd-Te = {VS/4.)acdTe = 2.87 A 
and take dun-Te = dcd-Te- The atomic structures are 
relaxed until the forces on each atom are smaller than 
0.02 eV/ A. 



III. NANOCRYSTALS WITH A SINGLE MN 
IMPURITY 

A. Energetics and Geometry for Mn Atom in 
Different Positions 

We study NCs doped with a unique Mn impurity sub- 
stituting a single Cd atom. The NC geometry is plotted 
in Fig. [T](a). The cation sites of the NC are distributed 
in three sets, and labeled as "I" for the sphere center, 
"11" for other inner positions and "III" for outside posi- 
tions. The substitutional energy is the energy difference 
AEi of the following reaction: 

CdTeNC + Mn+^ CdTeiMn NC + Cd+^ (1) 
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FIG. 1: (Color online) (a) Nanoparticle of (Cd,Mn)Te with 
107 atoms: Cd atoms with light gray, Te atoms with dark 
gray (green) and passivating H"^ atoms with small spheres. 
The cation sites are distributed in three groups; labelled as 
"I" for the center, "11" for inner atoms and "III" for nearest 
neighbors to surface atoms, (b) Substitutional energy A^i 
for a single Mn impurity. We are referring to the ions in 
solution. The higher is the energy the more stable is the Mn 
atom in such position. 



Due to the preparation of NCs in solution, we take as 
source systems for Mn and Cd their respective ions. The 
substitutional energies are shown in Fig. [1] (b) for the 
positions I, II, and III. All the energies are endothermic. 
This means that the synthesis of these NCs requires high 
temperatures as shown in experiments^^. The changes 
in AEi show opposite trends to the QD calculated total 
energies as a function of Mn site. The lowest-energy sites 
belong to group III with Mn atoms close to the surface; 
the outward sites are the most stable positions for the 
impurity. This outward stability is consistent with other 
results for magnetic impurities in other semiconductor 
nanoscale objects, such as nanowires^^. 

In the output geometry, the nearest-neighbor distances 
to the Mn impurity are about 2 % smaller than the Cd- 
Te distance for the undoped NC. This contraction effect 
depends strongly on the exchange-correlation approach. 



The neighbor Mn-Cd distances at the GGA level con- 
tract by more than 4 % respect to the bulk Cd-Te dis- 
tance. However, when relaxing the geometries with the 
-\-U scheme the distances expand closer to the unrelaxed 
input Cd-Te distances. This result could justify the use 
of unrelaxed Cd-Te distances in other works about Mn 
doping, but in principle this agreement seems fortuitous. 



Magnetism and Electronic Properties 



We analyze next the origin of magnetism as we are 
dealing with Mn, that typically is a magnetic element. 
The total magnetic moment associated with the QD is 
5 /iB, as expected since the Mn dopant introduces five 
spin- up electrons. The local magnetic moment in Mn 
impurity is nevertheless smaller, 4.65 /j^b- It differs sub- 
stantially from the bulk value, 4.21 /i^, for Mn in bulk 
CdTe as given in a recent work^^. The spin density is 
spatially plotted in Fig. [2l The negative areas of spin 
density in the Mn nearest-neighbor Te sites integrate to 
small magnetic moments about —0.01 fiB- Such Te atoms 
were non-magnetic before doping. Other radii for the in- 
tegration spheres different from the default Wigner-Seitz 
radii change only slightly these local magnetic moment 
values. Note that the -\-U approach affects the charge 
distribution and thus modify the local moment on Te, 
which changes in sign when improving the d-level de- 
scription. This behavior is in disagreement with other 
previous results concerning CdTe:Mn bulk^^, where the 
Te and Mn moments are found to be parallel. The Mn 
magnetic moment is lower than 5 fiB because of the sp- 
d hybridization between the localized 3d electrons of Mn 
impurity and the delocalized s- and p-type valence states 
of CdTe host. 

Before studying the sp-d hybridization, lets comment 
the gaps using the local density of states (DOS) as given 
in Fig. [3l Quantum confinement in QDs produces a blue 
shift of the gap and all the related optical properties, such 
as excitons. The gaps calculated using LDA or GGA ap- 
proximations are well known to underestimate the exper- 
imental values. The calculated gap of CdTe bulk is 0.69 
eV which is lower than in experiments, 1.48 eV. Our cal- 
culated HOMO-LUMO gap in the Cd-Te QD is 2.59 eV, 
which shows the predicted blue shift. All these findings 
are well established in semiconductor QDs. However, the 
HOMO-LUMO gap of the CdTe dot doped with Mn and 
within the -\-U approach remains almost constant to 2.62 
eV [see Fig. [3] (a)] when Mn is in the center. We note 
that the GGA approach understimates the gap and be- 
comes 2.04 eV. Moving Mn to other positions changes 
this gap value within a tenth of eV, which is negligible 
to be commented [see Fig. [3] (b)]. 
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FIG. 2: (Color online) (a) Spin-polarization isosurfaces when 
the Mn impurity is placed at the center for cuts at ±0.01 
jib/ A^. Positive values are light gray (red); negative val- 
ues, dark gray (blue), (b) Spin density in a plane defined by 
the central Mn impurity and two nearest-neighbor Te atoms. 
For the sake of clarity, spin-up density is chopped at 0.015 
/J.B A~^. The spin-down density region yields on the Te sites 
a magnetic moment of -0.01 /is, antiferromagnetically cou- 
pled to the Mn magnetic moment (4.65 /jb)- 



1. Origin of Mn QD-Exchange Coupling 

The sp-d hybridization is studied using the local den- 
sity of states (LDOS) projected into the orbitals of Mn 
impurity and of nearest-neighbor Te atoms. The main 
peak of Mn (i-states appears below the p part valence 
band. Due to the tetrahedral symmetry in CdTe lattice, 
the five spin-up Mn electrons are split into a triplet of 
t2 symmetry and a doublet of e symmetry. We are in- 
terested in the magnetic coupling between Mn and NC 
states equivalent to the band-edge states. The projected 
DOS around valence- and conduction-edge states are 
shown in Fig. [3l The conduction-edge states of the host 
compound show 5-type character and do not hybridize 
with the 3d orbitals of central Mn, or hybridize sligthly 
for off-center Mn. Hence, the s-d exchange interaction 
arises mainly from Coulomb repulsion and the Pauli ex- 
clusion principle^'^'^^, and originates the spin splitting 



of these states. This splitting is thus always ferromag- 
netic, and its exchange constant A^o<^ is positive. On the 
contrary, the valence band-edge states of the CdTe are 
p-type and allowed to hybridize fully with the Mn 3d or- 
bitals. Anyhow such p-d hybridization is small and yields 
an effective exchange mechanism of interaction between 
Mn atoms^^. This exchange is related to the opposite po- 
larization between the valence band-edge states and the 
Mn 3d orbitals. The p-d interaction originates the spin 
splitting of valence-band edge and it is always antifer- 
romagnetic, since all the 3d spin- up states are occupied 
and only jumps into the 3d spin-down states are avail- 
able. The corresponding exchange parameter NqP is thus 
always negative. 



2. Site Dependence of the Exchange Constants 

Now we compute the sp-d exchange interaction pa- 
rameters Noa and Nq^ following the expressions for 
bulks. They are defined in the standard mean- field the- 



ory as 
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The number A^o is the cations per unit volume^^. The 
differences AE^'^ = £^^'^(spin down) - £^^'^(spin up) are 
the spin splittings of the conduction (c)- and valence {v)- 
band edges. They can be extracted from the density of 
states depicted in Fig. [3l The x value is the fractional 
dopant concentration^^, and {Sz) = | is the average z 
component of Mn spins. 

Following eqn (2) we calculate Nool and A^o/^ as a func- 
tion of the impurity position. They are shown in Fig. IH 
When the Mn atom occupies the NC center, Noa = 0.41 
eV, No^ = -0.62 eV and No{a - f3) = 1.03 eV. The 
dopant concentration is x = 1/19 ~ 0.05 and these val- 
ues are comparable with the experimental ones obtained 
for bulk Cdo.95Mno.05Te: A^oQ^ = 0.22 eV, NoP = -0.88 
eV and No{a - f3) = lAO eV^^ For NCs, the absolute 
values of No/3 and Noa are close to each other. When 
the impurity is located off-center, the exchange constants 
are smaller than the bulk values and they tend to zero 
as the impurity approaches the surface. This finding is 
explained because the Mn bond is expected to be less 
covalent when Mn is located near the surface, where it 
is less bulk-like. This site dependence of the exchange 
constants seems to have implications for random distri- 
bution of NCs. When there is a collection of QDs, and 
each containing a single Mn impurity randomly placed, 
the average No{a — j3) would be significantly reduced 
in comparison with the bulk intrinsic Zeeman splitting. 
This decrease fits in the previous theoretical results^^^^^. 
The decrease of the exchange splitting is correlated with 
a smaller CdTe dot density around Mn. Thus, a further 
density analysis of band-edge states will be needed. 
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FIG. 3: (Color online) Projected DOS in Mn impurity, d- 
type, and in the nearst-neighbor atoms Cd and Te, s-type in 
black (black) and p-type in gray (orange). In panel (a) Mn is 
placed at the QD center; and panel (b), Mn at site II. Insets 
widen the conduction edge states. The dashed vertical lines 
separate the edges of valence band and conduction band. 



3. Densities of HOMO-LUMO States, in and off- Center 
Positions 

Though the detailed wavefunctions, as we have seen, 
are not required in the study of exchange interaction pa- 
rameters, we must clearly attribute them to the electron 
and the hole effective states which reflect the character of 
the conduction and valence states in the NCs doped with 
Mn. Our idea is to look at the High Occupied Molecular 
Orbital (HOMO) as a representation for the hole; and at 
the Lowest Unoccupied Molecular Orbitals (LUMO), for 
the electron. They are shown in Figs. [5]and[6l 
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FIG. 4: (Color online) Exchange splitting of band edges (left) 
and the Noa and Nq/S exchange constants (right) as a func- 
tion of the position (below) and the distance (above) between 
the impurity and the NC center. Positive values are LUMO 
splittings; and negative values, HOMO splittings. 



We do not include here the HOMO-LUMO states of the 
undoped NC because they are indistinguishable by simple 
eye inspection from the up wavefunction in Fig. [5] (a) and 
(c) for Mn in the dot center. The main contribution to 
the up HOMO comes from the nearest-neighbor Te atoms 
and is larger than from any other Te atom. Although 
more delocalized, the up LUMO has larger contributions 
in the Cd atoms and in the center site, which is Mn or 
Cd for doped or undoped dots respectively. These spatial 
distributions reflect the Te and Mn local DOS character 
for the HOMO and LUMO commented in the previous 
sections. 

However, the HOMO and LUMO down states are dif- 
ferent from the undoped NCs . For the down states we 
see that the Mn placed in the center expels charge, which 
is quantified by integrating it around the sphere center 
in Figs. [5] (e) and (f). This effect is much larger for the 
HOMO than for the LUMO, see insets. This difference 
can be explained partially because the HOMO state is 
occupied, and mainly because the LUMO state does not 
hybridize with the Mn states, as shown in the previous 
DOS plot [Fig. [3], which means that it has lower inter- 
action with the d-electrons of Mn. The down HOMO 
state undergoes the strongest change in comparison with 
its undoped counterpart because the undoped HOMO is 
mainly in the nearest-neighbor Te to the central Mn, with 
which the Mn must couple antiferromagnetically. 

To end this section we comment the density of the 
HOMO and LUMO states for off-center Mn. They are 
plotted for the position II in Fig. [6l There is larger mix- 
ing in the LUMOs which makes that both up and down 
states remove charge from the Mn neighborhood. The 
general form of the rest of the LUMOs resembles those 
of the sphere without Mn. The HOMOs with off-center 
Mn suffer larger differences with respect to the undoped 
case, specially for the down component. Anyhow, we see 
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FIG. 5: (Color online) Orbital densities in the Mn-centered 
sphere for the Highest Occupied Molecular Orbital (HOMO): 
(a) and (b) for up and down states. Also the densities for 
the Lowest Unoccupied Molecular Orbital are shown in (c) 
and (d) panels, up and down respectively. The density cut 
is around one third of the maximum value. The integrated 
charges around nanoparticle centers are given in (e) and (f) for 
HOMOs and LUMOs, respectively. Their up-down differences 
are plotted as dashed-dotted lines in gray. The Cd Wigner- 
Seitz Radius i?ws and the nearest-neighbor distance d nn are 
also shown with vertical dotted lines. 



a depletion of charge around central Mn mainly for the 
down part. These larger differences for the HOMO are 
due to two reasons: (i) the stronger mixing of states in 
a structure with lower symmetry for Mn off-center, and 
(ii) the small value for the undoped charge density out 
of the dot center. The latter means that when Mn is off- 
center, the states can suffer naturally stronger perturba- 
tions. Such is the case for the HOMO-LUMO states with 
off-center Mn. They show larger voids or lower values of 
the density around Mn. 



C. Analysis of Valence-Band Edge States and 
Occupied Mn Levels 

Exchange splitting of the conduction edge states has 
been discussed at some length in the previous sections. 
Here, we look again at the splitting of the valence-band 
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FIG. 6: (Color online) HOMO-LUMO densities for Mn off- 
center in position H. Panels (a)-(d) are organized as in the 
previous figure. The crossed atom denotes Mn. 



border by looking at the nearest border states. This anal- 
ysis is useful to clarify the breaking of symmetry in CdTe 
levels and Mn d levels due to the Mn displacement from 
the central position. 

The valence edge states for Mn in position I are shown 
with empty dashes in Fig. [71 (a) for up and down spins. 
The states are three-fold degenerate and we have com- 
mented at some length about their characteristics in the 
previous section. In this panel we show also the valence 
edge states for Mn in position II along the x direction 
of displacement. The off-center Mn breaks the degener- 
acy of these three states. We have identified the global 
characteristics of these states, although sometimes they 
have a strong mixed character. The Py and orbit als, 
perpendicular to Mn shift, have a smaller splitting. The 
global Px suffers the largest splitting as it is parallel to 
the Mn displacement direction. These electronic levels 
for Mn in position II explain naturally why the spin split- 
ting of up-down valence band edges is reduced to half 
value, respect to Mn at the center. 

These splittings are correlated to those of Mn levels, 
plotted in Fig. [3 (b). For central Mn the occupied d 
levels are grouped in two sets: the first has dx>'^_yii and 
d^ii^ the second gets the cross d levels x'y' ^ y' z' ^ and z'x' . 
These cross orbitals have lower energy which is typical for 
a Cd vacancy Vacd interacting with the tetrahedrally 
split Mn — d electrons [see branching scheme in Fig. [8]. 
The removal of a Cd atom creates a vacancy Vacd- The 
t2 levels of the Vacd hybridize with the ^2 states of Mn 
atom and produce the bonding t\ and antibonding ^2 
states. As the Mn states are very low in energy the 
states resembles much the Mn — 12 states and cross over 



the Mn — e states. While the ^ states remain close to 
the QD semiconductor states. For Mn in position II, the 
d levels shift upwards in energy, typical for the smaller 
CdTe density around Mn, as seen in the section about 
site dependence of exchange constants. The e levels re- 
main degenerate. The ^2 states {x^y\ y' z' ^ and z'x') split 
although they still are groupable together. The orbital 
z'x' orthogonal to the Mn displacement has the highest 
energy. The orbit als x'y' and y' z' remain at lower en- 
ergy, the lowest energy for the x'y' orbital with the lobes 
in the Mn displacement direction from the center. 

We want to stress finally that the Mn has a much 
larger exchange splitting of several eV that polarizes the 
CdTe levels in much smaller amount, in the order of sev- 
eral hundredths of eV. This spin polarization of CdTe in 
NCs depends on the Mn position. Although the splitting 
of CdTe levels that interact strongly with Mn is almost 
independent of position (P^c), the other levels {Py and 
Pz) remain almost unaffected. This finding means that 
the exchange splitting of the up-down CdTe valence edge 
amounts to nearly less than half of the value with central 
Mn. 



IV. NANOCRYSTALS WITH TWO MN 
IMPURITIES 

We investigate (Cd,Mn)Te NCs doped with two Mn 
impurities. The substitutional energy differences A£^2 
are calculated for the following reaction: 



Spin up 



CdTe NC + 2Mn+^ 



CdTeiMn NC + 2Cd 



+2 



(3) 



The calculated differences against the Mn positions 
and the distance between Mn atoms are displayed in Fig. 
m When both Mn impurities occupy sites I-II and II-II, 
the energies A£^2 are larger than for positions I-III and 
II-III. This difference is about 0.5 eV. The substitutional 
energy of Cd decreases with the presence of Mn atoms 
at surface sites III. This trend follows the previous en- 
ergy differences /S.Ei for a single Mn in Fig. (TJ where 
position III is the most stable. The interaction energy is 
the difference between the previous A£^2 Mn-Mn energy 
and the sum of single Mn energies A^i in the NC. We 
plot also the interaction energy (A) between the two Mn 
impurities in Fig. [9](b). The changes of the interaction 
energy A due to Mn positions are an order of magnitude 
lower than those of the substitutional energy A£^2- The 
interaction A is 32 meV for Mn-Mn atoms in sites I-II 
and 21 meV in III-III. It seems that the closer the Mn 
impurities the larger their interaction energy. However, 
the dependence of A with the distance between the two 
Mn impurities is almost flat except when Mn atoms are 
next-nearest neighbors or close to the surface. 

To a lesser extent these energies A£^2 and A also de- 
pend on the magnetic configuration between Mn mag- 
netic moments. We consider two magnetic coupling cases 
for Mn moments: ferromagnetic in which fiQD = 10 /i^, 




• ft 

x'y', y'z\ z'x 



FIG. 7: (Color online) (a) Valence states close to the Fermi 
energy for Mn in center with empty lines; and for Mn off- 
center in position II with full lines. The global densities are 
given nearest to the corresponding states. The notation con- 
cerning atom labeling and density cuts follows the one given 
in previous figures. The geometrical insets show the reference 
systems. The degeneracy gets lower for the off-center posi- 
tion, however the quantum dot states Py and Pz are nearly 
degenerate, (b) Mn d states for center and off-center posi- 
tion. We see the splitting between the triplet t2 and doublet 
e states [see next figure for their order explanation]. 
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FIG. 8: (Color online) Branching diagram for a substitutional 
Mn atom in CdTe dot. The states of substitutional Mn atom 
(center) are originated in the hybridization of vacant Cd states 
Voicd — t2 (right) with Mn-d states (left). 



and antiferromagnetic in which jhqd = /i^- Indepen- 
dently of the Mn positions within the QD, the Mn mag- 
netic moments in the lowest-energy state are coupled an- 
tiferromagnetically. This is shown in Fig. [10] using the 
d-d exchange constant for various Mn-Mn positions in- 
stead of the FM-AFM energy differences. In fact, we 
calculate the d-d exchange constant from the FM-AFM 
energy differences for various Mn-Mn positions. These 
constants J^^ are depicted in Fig. [10] (a). The largest 
value in modulus is roughly 1.5 K for Mn atoms on sites 

I and II. It is about four times smaller than the abso- 
lute value of the experimental exchange constant in bulk 
Cdo.95Mno.05Te, 6.1 K. In addition we can see the depen- 
dence with the Mn-Mn sites. For Mn-Mn intermediate 
distances, the exchange constants decrease even further 
and they are following the AFM-FM splitting. Note that 
for the Mn atoms in the center and surface positions this 
splitting seems to be very small. 

We have also examined the (i-levels for several of these 
Mn-Mn configurations, although they are not given here. 
Due to symmetry in QDs, the d-levels in position I and 

II are close, but in different energy panels [as seen in 
previous Fig. [7]. This means that the symmetry con- 
siderations are more relevant in QDs to understand the 
Mn-Mn interaction that higher-order mechanisms, such 
as superexchange or double exchange. We must remem- 
ber that such mechanisms are important for DMS bulks 
where the levels are at the same energies, due the equal 
density around Mn atoms. These symmetry considera- 
tions and the smaller dot density around Mn atoms ex- 
plain the lower coupling J^^ constants. 

The A^o<^ and No/3 exchange constants are shown in 
Fig. [To] (b). They are calculated for the ferromag- 
netic (circles) and antiferromagnetic (squares) states for 
several Mn-Mn positions; the dopant concentration is 
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FIG. 9: (Color online) (a) Substitutional energy of two Mn 
positions, (b) Interaction energy A for various Mn-Mn pairs. 
It shows the difference between the previous total Mn-Mn 
energy and the sum of single Mn energies in the NC. 



X = 2/19 ~ 0.1. As in the previous section, the exchange 
constants get closer to the bulk experimental values when 
one of the Mn impurities occupies the NC center (I-II and 
I-III), and they tend to zero when the Mn ions separate 
from each other and approach the surface. The exchange 
constants for antiferromagnetic cases with high symme- 
try are also zero because the down and the up electrons 
are fully degenerated. When comparing Fig. [4] with Fig. 
[TOl we conclude that the largest exchange values are ob- 
tained for NCs doped only with a Mn impurity in the 
central site. The relative orientation of Noa and NqP is 
the same as before. However their values decrease by at 
least 0.2 eV. This difference can be ascribed to the sec- 
ond Mn atom that cancels partially the coupling to the 
conduction and valence bands in the neighborhood of the 
first Mn. 



A. Valence Edge States 

The valence states around the edges allow us, like in 
the case of a single Mn impurity analysis, to understand 
the origin of the exchange in the CdTe states. As it has 
been done before, they are plotted for the asymmetric 
Mn-Mn configuration I-II in Fig. [11] and for the sym- 
metric one II-II in Fig. [12] However, due to magnetic 
coupling there are more options than in one Mn impurity 
case. Thus the panels (a) denote the FM Mn alignments; 
and the panels (b), the AFMs. 

We shall see that for the results of both figures, a per- 
turbation theory of the central Mn case rationalizes also 
the level breaking as due to symmetry and magnetic con- 
figuration. We therefore start the discussion by consid- 
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FIG. 10: (Color online) The (a) d-d exchange constant to- 
gether with the (b) sp-d exchange constants (A^o^ and Nol3) 
for several Mn-Mn positions. The J^^ are given in K to com- 
pare with the bibliography. Exchange splittings follow the 
notation of Fig. 4: Circle marks with dotted lines (green) 
denote FM alignments; and square marks with dashed lines 
(blue), AFM alignments. 



ering the AFM case I-II in Fig. [TT] (b). 

We see that the exchange splitting of the CdTe aligns 
with the Mn in position 1. Let us consider that a Mn in 
position II perturbs the levels to first order, as in previous 
one Mn impurity. Now, we can focus on the change in 
splitting of the three-fold degenerate levels using the pre- 
vious results. In this way, the down CdTe valence-band 
edge levels are Px at low energy, and nearly degenerate 
Py and Pz at higher energies. For the up levels this order 
is reversed. This order is also followed by our computa- 
tions. The expected and calculated levels are in excellent 
agreement and thus confirm that a perturbative approach 
applies to the splitting of CdTe levels by central Mn. 

However, the case AFM I-II we have discussed, in 
which the perturbative approach can be applied, is not 
the only case of physical interest. In the FM case. Fig. 
[TTI (a) for example, both Mn reinforce their splitting of 
the CdTe states. The value is close to the sum of both I 



and II splittings for the Mn impurities on their own. We 
shall not go into further detail here. 

For the case II-II we only draw attention to the fact 
that there is strong interaction for one up Mn impurity 
with the CdTe down states. We can see the empty space 
of up (down) orbitals around the down (up) Mn in Fig. 
[T2I (b) as it expels the surrounding charge. We want to 
stress here that for the FM case in Fig. [12] (a) we need to 
go beyond the previous simple perturbative analysis. In 
spite of the strong perturbation, it means that it is more 
difficult to interpret the data as coming from a perturba- 
tion due to a second impurity. It remains true that there 
is a clear increase of the splitting of Px orbital aligned 
along the corresponding Mn-Mn direction. 



V. CONCLUSIONS 

In summary, we have investigated the electronic and 
magnetic properties of (Cd,Mn)Te nanocrystals of spher- 
ical shape with the density functional theory. The em- 
bedded Mn impurities substitute Cd atoms in the zinc 
blende lattice. The QDs are ~ 2 nm in diameter, cen- 
tered on a cation site (Cd), and they have a total of 107 
atoms. The Cd and Te dangling bonds on the surface 
are passivated by pseudohydrogen atoms (H*), which are 
the simplest model for a real organic ligand. We have 
studied two doping situations: NCs including one and 
two Mn impurities. 

In the first case we find that the configuration with 
the Mn atom near the surface presents the lowest energy. 
The Mn impurity introduces five c?-type spin-up electrons 
within the NC, thus the total magnetic moment associ- 
ated with the QD is 5 jub^ as expected. The local Mn 
magnetic moment is 4.65 fiB- It is smaller than 5 /i^ 
because of the s,p-d hybridization. Also the Mn- nearest- 
neighbor Te sites show small magnetic moments antifer- 
romagnetically coupled to the Mn moment. Furthermore, 
we look at the hybridization in the density of states. We 
calculate the sp-d exchange constants for various Mn lo- 
cations. The exchange constants are comparable to the 
bulk ones when the impurity is at the center, and they go 
to zero when the Mn is close to the NC surface. Then, 
for central Mn we introduce spatially the HOMO and 
LUMO states to verify that the largest change is for the 
down HOMO, which also justifies the antiferromagnetic 
alignment of the nearest-neighbor Te atoms. Addition- 
ally, we check the role of symmetry by looking at the 
CdTe levels near the valence-band edge. This allows us 
a clear interpretation of the origin of exchange splitting. 
The Px levels in the Mn diplacement direction suffer the 
stronger hybridizaton, while the Py and Pz levels remain 
almost unsplit. Thus, the exchange contants are roughly 
half of the central Mn case. 

In the second case we calculate the total ground-state 
energy as a function of Mn-Mn positions and the mag- 
netic configuration of their local magnetic moments. In 
the minimum-energy state the Mn dopants are on the 
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FIG. 11: (Color online) Valence states closest to the Fermi 
energy for Mn in configuration I-II with full lines for (a) FM 
and for (b) AFM couplings. The empty dashes refer to a pre- 
vious case with Mn in a central position. The global densities 
are given nearest the corresponding states. The notation con- 
cerning atom labeling and density cuts follows the one given 
in previous figures. At the bottom, a scheme of Mn magnetic 
moments is shown. For the ground AFM case the valence 
states split according to center Mn in position I, however the 
degeneracy follows closely the off-center Mn in position II. 



FIG. 12: (Color online) Valence states closest to the Fermi en- 
ergy for Mn in symmetric configuration II-II with full dashes 
for (a) FM and for (b) AFM coupUngs. They follow the nota- 
tion of the previous figure. For off-center positions the valence 
states break the degeneracy even further. 



surface, as far apart as possible, and antiferromagneti- 
cally coupled. We calculate the sp-d and d-d exchange 
constants for various Mn-Mn locations. We find that the 
exchange constants A^o<^ and NqP are comparable to the 
corresponding bulk values when impurities occupy cen- 
tral bulk-like positions (I-II). They tend to zero as they 



separate from each other or approach the surface. As for 
the exchange |J^^|, its largest value is for the Mn ions 
placed in central bulk-like sites (I-II), and is four times 
smaller than in the bulk with similar Mn concentration. 
A similar analysis of states nearest the valence-band edge 
demonstrates that a perturbative approach seems to be 
valid, with exception of the II-II FM configuration. 
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